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NOTATION 


The notation used in the following sections is defined with the corre- 
sponding FORTRAN name used in the program indicated parenthetically. Figure 1 
shows the nomenclature used for the various inlet and nozzle locations. The 
values below are defined per engine and the drag coefficients are based on 
inlet capture area unless noted. The starred (*) items are required program 
inputs which are either user input or are supplied by another subroutine in 
the aircraft synthesis program. 


Symbol 

Code 


A 


area, m^, ft^ 

\uX^*ENG 

(AUAENG) * 

auxiliary systems area ratio 


(ABLEAC) 

bleed mass flow ratio 


(ABYPAC) 

bypass mass flow ratio 

A 

c 

(AC) 

inlet capture area (per engine), m^, ft^ 

^cc 

(ACC) 

area of exit nozzle (joint point between engine 
and fuselage) 



area of exit, m^, ft^ 

^EF 

(AEF)* 

engine face flow area (per engine), m^, ft^ 

^NG 

(AENG)* 

engine face total area (per engine), m^, ft^ 

^IT 

(AEXIT) 

nozzle exit area (per engine), m^, ft^ 


(ANOZT) 

nozzle throat area (per engine), m^, ft^ 

A 

o 

(AO) 

area of free-stream stream tube (per engine), m^, ft^ 

A /A 
0 c 

(AOAC) 

mass flow ratio of inlet (per engine), m^, ft^ 

A 

s 

(AS) 

projected frontal area of compression surface, 
m^, ft^ 


(AT) 

inlet throat area (per engine), m^, ft^ 


(ATD) 

inlet throat area (per engine) at 

*VENT^''c 

(AVEACD) 

ratio of engine ventilation flow area to inlet 
capture area (per engine) 


iii 



\eDGE^^c 

A 

y 


AD 


"AUX 


BL 




BP 


BT 


DIV 


INF 


DIV 


Cg or C. 


SP 


cc 

^ENG 

g 

D 

g 

h 

IPR 

L 


(AWAENG) * 
(AY) 


(CDAD) 

(CDAUX) 

(CDBE) 

(CDBP) 

(CDBT) 

(CDDIV) 

(CDI) 


(CPCS) 

(CS or CDADS) 


boundary-layer diverter area ratio 

projected frontal area of compression surface 
forward of point of normal shock impingement, 
m2, ft2 

drag coefficient 

supersonic spill additive drag coefficient 
auxiliary systems drag coefficient 
bleed drag coefficient 
bypass drag coefficient 
nozzle boattail drag coefficient 
boundary-layer diverter drag coefficient 
nozzle interference drag coefficient 
boattail drag coefficient based on 
pressure coefficient on diverter surface 
pressure coefficient on compression surface 
subsonic spill additive drag coefficient 


thrust coefficient 

(DCC) nozzle diameter at customer connect, m, ft 

(DENG) engine face diameter, m, ft 

acceleration of gravity, m/sec^, ft/sec^ 

(DEXIT) nozzle exit diameter, m, ft 

altitude, m, ft 

(IPR)* inlet pressure recovery code 

distance between normal shock position and inlet 
lip 


iv 




(XLVD) 

distance between normal shock position and inlet lip 


(XLNOZ) 

ratioed to inlet capture diameter 

^NOZ 

nozzle length, m, ft 

M 


Mach number 

in 


mass flow, kg/sec, Ib/sec 

“aux 


auxiliary systems mass flow, kg/sec, Ib/sec 



bypass mass flow, kg/sec, Ib/sec 

M 

cone 

(XMCONE) 

compression surface Mach number 


(XMDES)* 

inlet design Mach number 



exit Mach number 



exit mass flow, kg/sec, Ib/sec 


(XMEF)* 

engine face Mach number 

^XIT 

(XMEX) 

nozzle exit Mach number 


(XMT)* 

inlet throat Mach number 

M 

00 

(XMO)* 

free-stream Mach number 

^ENG 

(EN)* 

number of engines 

NPR 

(NPR) * 

nozzle pressure ratio 

P 


static pressure, N/m^, Ib/ft^ 

^E 


exit static pressure, N/m^, Ib/ft^ 

^Ses 

(PRDES) 

supersonic diffuser pressure recovery at 

**^SUB 

(PRSUB) 

subsonic diffuser pressure recovery 

^^SUP 

(PR) 

supersonic diffuser pressure recovery 


(PRTOT) 

total pressure recovery to engine face 

PSPIN 


cone surface pressure ratio 



total pressure, N/m^, Ib/ft^ 


(PTBLE) 

bleed exit total pressure, N/m^, Ib/ft^ 

Bleed 




V 



Bypass 

(PTBYP) 

bypass exit total pressure, N/m^, Ib/ft^ 



total pressure at engine face, N/m^, Ib/ft^ 

^TH 


cone static pressure at the throat, N/m^, Ib/ft^ 

^NOZ 

(PTNOZ)* 

nozzle exit total pressure, N/m^, Ib/ft^ 

P 

^TH 


total pressure at inlet face, N/m^, Ib/ft^ 

P 

t 

00 

(PTO)* 

free-stream total pressure, N/m^, Ib/ft^ 

p 

oo 

(PINF)* 

free-stream static pressure, N/m^, Ib/ft^ 

Q or q 
^00 

(Q)* 

free-stream dynamic pressure, N/m^, Ib/ft^ 

SFC 


specific fuel consumption, kg/N-hr, Ib/lb-hr 

S/D 

g 

(SODG)* 

nozzle spacing ratio 

^ef 

(SWING)* 

wing reference area, m^, ft^ 

T 


thrust, N, lb 

T 

S 

(FIP)* 

gross thrust per engine, N, lb 

■'t 


total temperature, K, R 

NOZ 

(TTNOZ)* 

nozzle exit total temperature, K, R 

8 

4^ 

H 

(TTO)* 

free~stream total temperature, K, R 



exit velocity, m/sec, ft/sec 



free-stream velocity, m/sec, ft/sec 

w 

a 

(WA)* 

engine airflow, kg/sec, Ib/sec 


(XCOYC) 

distance from cone tip to inlet face ratioed to inlet 
capture diameter 


(YC) 

inlet capture diameter, m, ft 



diameter of inlet centerbody at inlet throat, m, ft 

e 

(BETA) 

nozzle boattall angle , deg 


vl 



(DELPR)* 


(GAMMA) 

(LAMBDA) 

(RHO) 

(THETA) 

(THDIV) 


incremental pressure recovery correction 
isentropic constant 

angle at inlet lip between average direction of flow 
and longitudinal axis of inlet 

free-stream static density, kg/m^, Ib/ft^ 
cone half angle, deg 

boundary-layer diverter wedge angle, deg 

exit angle, deg (COSDE is cosine of exit angle in 
program) 



A SIMPLIFIED ANALYSIS OF PROPULSION INSTALLATION 


LOSSES FOR COMPUTERIZED AIRCRAFT DESIGN 
Shelby J. Morris, Jr., Walter P. Nelms, Jr., and Rodney 0. Bailey 

Ames Research Center 


SUMMARY 


A simplified method is presented for computing the installation losses of 
aircraft gas-turbine propulsion systems. The method has been programmed for 
use in computer-aided conceptual aircraft design studies that cover a broad 
range of Mach numbers and altitudes. The items computed are: inlet size, 

pressure recovery, additive drag, subsonic spillage drag, bleed and bypass 
drags, auxiliary air systems drag, boundary-layer diverter drag, nozzle boat- 
tail drag, and the interference drag on the region adjacent to multiple nozzle 
installations. The methods for computing each of these installation effects are 
described and computer codes for the calculation of these effects are furnished. 
The results of these methods are compared with selected data for the F-5A and 
other aircraft. The computer program can be used with uninstalled engine 
performance information which is currently supplied by a cycle analysis pro- 
gram. The program, including comments, is about 600 FORTRAN statements long, 
and uses both theoretical and empirical techniques. 


INTRODUCTION 


The design of advanced aircraft systems requires the consideration of many 
different tradeoffs and parameters to arrive at an optimum design for a 
particular requirement or group of requirements. One is the effect of inter- 
action between the aerodynamics and the propulsion of these systems. Propul- 
sion installation effects on high-speed aircraft can amount to 10 percent or 
more of the aircraft drag and can also degrade the propulsion thrust via inlet 
total-pressure recovery penalties and nozzle-flow penalties. These effects 
are significant in high-speed aircraft design, and thus require attention, 
even in early design studies. 

Tradeoff studies are usually done manually or, more recently, by many 
large computer programs with manual communication between them. As computer 
capabilities have increased, it has become possible to communicate between 
these disciplines within the computer in an automated or integrated fashion. 
This integration allows computation of the trajectory of the aircraft over its 
entire mission, thereby providing the ability to determine the effects of 
various parameters and to optimize the aircraft for specific requirements 
subject to various constraints. The method and computer code presented in this 
report is intended to supply the propulsion Installation losses as required 
in this process. The code is designed to work as part of a propulsion module 


in the framework of the Aircraft Synthesis Program, ACSYNT (fig* 2), which 
has been developed at the Ames Research Center (ref. 1). 


The purpose of this report is to document the methods and the computer 
code for propulsion installation losses as presently employed in ACSYNT. 
Limited example comparisons of calculations with data are made and areas of 
further research identified. It should be emphasized that, at present, the 
methods are preliminary in nature and further work is needed to improve the 
techniques and to perform additional correlations with data. 


PROGRAM PHILOSOPHY 


The purpose of the Propulsion Jnstallation (7alculation (PRINC) module is 
to compute the air induction system and nozzle/afterbody effects in the 
ACSYNT program. The procedures employed in the present subroutine are general, 
since the methods must be applicable to a variety of inlet, engine, and nozzle 
types over a broad range of Mach numbers and altitudes. An additional impor- 
tant requirement is that the calculations be very rapid, since installation 
losses are computed many times (over 1000) in a run of the ACSYNT program. 

Figure 3 shows a block diagram of the method. A modular approach is used 
so that future additions and improvements can be easily incorporated. Items 
computed include (1) inlet pressure recovery, (2) inlet size, (3) additive and 
spillage drags, (4) bleed and bypass drags, (5) auxiliary system drag, 

(6) boundary-layer diverter drag, (7) nozzle boattail drag, and (8) nozzle 
interference drag. In figure 3, those parameters listed inside the boxes 
are output from the various modules and those parameters listed beneath each 
box are required inputs to each module. 

There are varied accounting approaches for the aerodynamic propulsion 
system and propulsion system/airframe interaction losses. The method employed 
in the PRINC module is to charge all losses (listed above) to the engine thrust 
and specific fuel consumption (SFC) as indicated in figure 4, However, the 
total propulsion installation drag as well as the individual propulsion- 
related drags are computed separately so that any desired accounting method 
may be adopted by the user. An available option in the program is a multiply- 
ing factor for any or all of the propulsion Installation losses to adjust the 
level of these penalties at the user's discretion. 


DESCRIPTION OF METHODS 


This section documents the methods used in the propulsion Installation 
loss module (PRINC) and diagrammed in figure 3. It is assumed, for the inlet 
drag calculations, that the inlet is an axlsymmetric, external compression 
design and, for the additive drag calculation, that the surface pressures are 
for a cone of an average half angle of 20°. The drag coefficients computed 
in the following development are based on inlet capture area, except where 
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noted* The equations, derivations, and programming details are presented in 
appendix A. A FORTRAN listing of all the modules is included in appendix B. 

Inlet Pressure Recovery — The inlet pressure recovery is divided into two 
parts, the pressure recovery in the region ahead of the inlet face and the 
pressure recovery in the subsonic diffuser after the inlet face. The pressure 
recovery in the region ahead of the inlet face is estimated by the use of 
the standard AIA or Military Specification 5008B methods or by the assumption 
of normal shock pressure recovery (appendix A) . The pressure recovery versus 
Mach number computed by these three methods is shown in figure 5. 

The subsonic diffuser pressure recovery is estimated by the empirical 
method of Ball (ref. 2), which gives this pressure recovery as a function of 
the throat Mach number, the inlet lip bluntness, and the free-stream Mach 
number. For the present study, the inlet lip has been assumed to be sharp and, 
thus, the inlet subsonic diffuser pressure recovery is independent of lip 
bluntness or free-stream Mach number. Also, the geometric inlet throat Mach 
number is equal to the effective inlet throat Mach number as described in 
reference 2. 

A fourth method available in the program is to input the inlet total 
pressure recovery as a function of free-stream Mach number in tabularized 
form. 


Inlet Sizing -* The inlet face flow area is determined by a mass balance 
(conservation of mass) between the inlet face and the engine face. The mass 
flow at the engine face is determined by the requirements of the engine. The 
inlet face flow area is increased over that of the engine to allow for bypass, 
bleed, and powerplant ventilation mass-flow requirements. The free-stream 
stream-tube cross-sectional area is determined by a mass balance between the 
free stream and the inlet face. The inlet design Mach number is used to define 
the inlet capture area, which is equal to the free-stream stream-tube cross- 
sectional area at the engine’s maximum power setting. The inlet capture area 
is held constant at off-design conditions; however, the centerbody is allowed 
to move so that the inlet throat Mach number is held at some specified value. 

No check is made on the mechanical difficulty of achieving this variation. 

The key assumption in this analysis is that the inlet throat Mach number is 
constant. The programming details of this subroutine are included in appen- 
dix A. 

Additive Drag — The engine thrust is referenced to free-stream conditions. 
The loss in momentum of the airflow ahead of the inlet system must be accounted 
for in the bookkeeping system. This loss in momentum ahead of the inlet face 
is called "additive drag" and is a function of the inlet geometry, the free- 
stream Mach number, and the mass flow of the engine. 

The inlet additive drag is computed by a momentum balance between the 
inlet face and the free stream. The cosine of average flow angle (with respect 
to the inlet centerline) at the inlet face is assumed to be 1.0. The inlet is 
assumed to be external compression (that is the normal shock is outside of the 
cowl lip). The inlet throat Mach number is held constant at some specified 
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value . 
can be 


The inlet geometry is assumed to be axisymmetric. The additive drag 
computed (ref, 3) from 


C 


AD 



A P 
^ c « 


^TH F 


TH 



+ 1 ) 


cos X 


+ 


Ar “ A' 


TH 


cone 


- 1.0 - 



+ C 


S 


The cone pressure calculation uses a polynomial approximation presented 
by Lighthill (ref, 4). The subsonic spillage effect Cg is computed using 
an empirical technique described by Sibulkin (ref. 3). A complete description 
of the method is included in appendix A. 


Bypass Drag — In high-Mach-number aircraft design the inlet is usually 
sized at the maximum design Mach number. During off-design operation at lower 
Mach numbers, the inlet usually has the capacity to supply an excess airflow 
to the engine. This excess airflow must be either taken onboard the aircraft 
and passed (bypassed) around the engine or diverted (spilled) around the inlet 
system. 


The bypass drag is computed from a momentum balance between the free- 
stream and the bypass exit. The bypass exit nozzle can be either sonic or 
fully expanded. After considerable simplification (see appendix A), the 
momentum balance yields 


7aT ' 2 ^ “ 

c ^ 




“e 

““ ®E 5” 


\1 + 0.2m|/ J 


(cos M„ 

/I + 0.2Me2\^ I 

1 + 

(0.7M^ 

' 00 

\1 + 0.2M^ / 1 

00 

_(P^ /P^ ) Vl + 0.2M^2j J) 

E 00 


where y is assumed to be 1.4. If It is asstimed that the bypass exit nozzle 
is sonic, then 

^ = 1.0 

If it is assumed that the bypass exit nozzle is fully expanded, then 

Pe “ 


- [5(P^ /pjO.286 _ 
E 
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The bypass exit pressure recovery is assumed to be a fraction of the inlet 
total pressure recovery (to the engine face) . Typical values for this fraction 
are 


E o 


KP /P 

EF ^ 


where 0.3 < K < 0, 7. 

Bleed Drag — The inlet compression ramp or cone for typical supersonic 
inlet designs often have a considerable length exposed to an adverse pressure 
gradient. This can create a boundary layer which is thick enough to cause 
losses in engine performance. The problem is particularly acute in regions 
where a shock wave interacts with this boundary layer. In order to maintain 
efficient engine performance, part of the boundary layer is removed on these 
compression surfaces in some inlets, and it is necessary to account for the 
momentum loss of this bleed flow. A momentum balance between the free stream 
and the bleed exit yields an expression similar to the bypass drag formulation. 
The bleed exit can be assumed to be either sonic or fully expanded. The 
momentum balance yields 



where y is assumed to be equal to 1.4. If it is assumed that the bleed exit 
nozzle is sonic, then 


Me = 1.0 

If it is assumed that the bleed exit nozzle is fully expanded, then 


= [5(P /PJO'286 _ i]0.5 

E 

The bleed exit pressure recovery is assumed to be a fraction of the inlet 
total pressure recovery (to the engine face) . Typical values for this fraction 
are 




KP /P 
EF 


t 

OO 


where 0.3 < K < 0.7. 

A complete derivation of these equations is contained in appendix A. 
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Auxiliary Systems Drag — The auxiliary systems drag accounts for the air- 
flow taken into the aircraft for systems cooling and auxiliary power generation. 
Many aircraft have small auxiliary inlets mounted at some convenient place to 
serve this purpose, and the drag created can be significant. It is assumed 
that the total momentum of the flow into these systems is lost. Therefore the 
auxiliary system drag is 


AUX 


”^AUX^“ 


P A 


V" 


” AUX 
4 P V^A 

2 00 00 Q 


AUX 


where 


A^Ux^^c ratio of the auxiliary system inlet capture area to 


aircraft inlet capture area. Typical values for this quantity range from 
0.005 to 0.01. 


Boundary-Layer Diverter Drag - In many inlet installation systems, the 
inlets are located close to the aircraft’s larger components (i.e., wings, 
fuselage) which generate regions of low momentum ahead of the inlet. The 
ingestion of these boundary layers into the inlet creates a nonuniform flow 
distribution which can cause considerable performance degradation in the 
engine. This problem has been avoided by the addition of a ramp (a plow) 
between the inlet and the boundary-layer generating surface. The turning of 
the flow in these systems adds drag to the aircraft, which must be accounted 
for. A fit of data (refs. 5 and 6) yields 


c = ^ ^DGE 

20 A^ 


; M > 1,55 

oo — 


= 0.499 ; 0.95 < M <1.55 

20 A - oo - 


0.499 


^oo " \fEDGE 


(0.95- 0.80) 20 A 


;• 0.80 < M <0.95 


— OO — 


=0.0 ; M < 0.8 

Details on the data and a comparison with the fit are given in appendix A. 

Boattail Drag — The boattail drag on the airframe back to the point where 
the nacelle and engine are joined (see fig. lb) is calculated as part of the 
aircraft drag. The boattail drag on the portion of the engine which includes 
the engine nozzle after this joint is charged to the engine performance in 
the present accounting system. The boattail drag estimation method used is an 
empirical technique developed by Ball (ref. 2) from wind-tunnel data on iso- 
lated boattail nozzles. The nozzle interference drag described in the next 
section corrects this for installations of more than one engine. The boattail 
drag is based on the area at the point where the engine is joined to the 
airframe. The formulation is for an engine nozzle pressure ratio (engine exit 
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total pressure to free-stream static pressure) of 2.5; however, correction 
terms are included for different nozzle pressure ratios. The engine nozzle 
exit area is computed from the engine thermodynamic data. The boattail angle 
is computed from the engine diameter and the assumption that the length of the 
boattail is equal to the engine diameter. It is also assumed that the diam- 
eter of the boattail at the connection point between the engine and aft 
fuselage or nacelle is 10 percent greater than the engine diameter. A complete 
description of this procedure is included in appendix A. 

Nozzle Interference Drag — The nozzle interference drag accounts for the 
drag on the base area between multiple nozzles. The independent variables are 
free-stream Mach number and nozzle spacing ratio S/Dg (ratio of the distance 
between nozzle centerlines to nozzle exit diameter). The calculation technique, 
developed by Ball (ref. 2) from wind-tunnel data, estimates the ratio of the 
drag due to nozzle interference divided by ideal gross thrust at a nozzle 
pressure ratio of 2.5. This value is corrected to a drag coefficient based on 
inlet capture area. A complete description of this computation is included 
in appendix A. 


EXAMPLE CALCULATIONS 


This section presents example computations from the PRINC module of 
typical installation drags, net propulsive thrust, and specific fuel consump- 
tion values. After PRINC module calculations of inlet mass flow and propulsion 
installation drags for a simulated F-5A are presented, these results are then 
used to determine the overall installed thrust and SFC of an ACSYNT simulated 
F-5A. Comparisons are made of these results with F-5A flight test data. 


Mass Flow Summary 

The effect of Mach number on engine mass flow ratio A^/A^ for the PRINC 
module simulated F-5A is presented in figure 6. Note that the F-5A has no 
bleed or bypass. The spillage mass flow is the difference between Aq/A^ = 1.0 
and the Aq/A^ set by the engine (plotted). This difference would be much 
larger for an aircraft with a higher inlet design Mach number Mqes- Th® 
method is capable of handling bleed and bypass in the manner described in the 
section on bypass and bleed drag. 


Total Installation Drag 

Figure 7 is an example PRINC module calculation of the installation drag 
coefficients based on wing reference area as a function of Moo for a 
simulated F-5A inlet system. The total installation drag coefficient is shown, 
as well as the various components for maximum afterburning (A/B) and military 
power settings. For this same inlet system, the effects of engine throttling 
at M - 0.9 and 1.2 are shown in figures 8a and b. 
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Net Propulsive Thrust Correlation 


A comparison of the thrust calculated by the ACSYNT propulsion subroutine 
and the PRINC module with data determined from F-5A flight tests is shown in 
figures 9a and b for maximum A/B and military power settings. The results are 
presented for two engines over a range of Mach numbers at 10 973 m (36 000 ft). 
The upper portion of each figure compares the uninstalled thrust from the 
ACSYNT propulsion module with corresponding values from the J-85-GE-13 engine 
specifications (ref. 7). Both thrust values are based on the AIA standard ram 
recovery schedule. The table shows the percentage difference between the 
calculated results and data for selected Mach numbers; that is, 


Calculated-Actual _ 

X XUU 

Actual 

The lower portion of the figure shows a comparison between the installed 
thrust calculated by the ACSYNT propulsion subroutine with corrections calcu- 
lated by the PRINC module and flight-test modified data from reference 8. The 
PRINC module calculations include corrections for a pressure recovery schedule 
based on a corrected airflow of 20.4 kg/sec (45 Ib/sec) (ref. 9) and for the 
following installation losses — additive drag, auxiliary systems drag, boundary- 
layer diverter drag, and nozzle boattail and interference drags. Bleed and 
bypass drags are zero. Exactly what corrections are included in the flight- 
test modified data of reference 8 is not clear, but it is suspected that losses 
for the boundary-layer diverter and the nozzle are not included. This would 
account for some of the overcorrection by the PRINC module. With a few 
exceptions, the percentage differences for both power settings are within 
10 percent. 


SFC Correlation 

Figures 10a and b show comparisons between specific fuel consumption 
values from the ACSYNT propulsion subroutine and the PRINC module and data 
determined from F-5A flight tests. These comparisons correspond to the thrust 
correlations shovm in figures 9a and b. As with thrust, the percentage 
differences are generally within 10 percent. It should be noted that the F-5A 
flight-test evaluation may use a different method of bookkeeping, which could 
account for some of the differences. 


CONCLUDING REMARKS 


A simplified method has been presented for computing the installation 
losses of aircraft gas-turbine propulsion systems. The program employs rapid 
and sufficiently accurate estimating procedures suitable for use in computer- 
aided conceptual design studies of aircraft systems over a broad range of Mach 
numbers and altitudes. The items which can be computed are: inlet size and 

pressure recovery, additive drag, subsonic spillage drag, bleed and bypass 
drag, auxiliary air systems drag, boundary-layer diverter drag, nozzle boattail 
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drag, and the interference drag on the region adjacent to multiple nozzle 
installations. The methods for computing each of these installation effects 
have been described and compared with either data or the results of more 
elaborate computing procedures. Finally, a comparison of the overall results 
of the method with F-5A performance specifications indicates an accuracy within 
a out 10 percent in installed thrust and specific fuel consumption. This is 
considered sufficiently accurate for computerized design at the early stages 
of vehicle definition. 
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APPENDIX A 


DEVELOPMENT OF PROGRAMMED EQUATIONS 


This appendix contains a brief development and description of the 
equations that are used in the PRINC program- The equations are presented by 
subroutine. 


INLET PRESSURE RECOVERY 
(MODULES PRSUBS AND PRINL) 


This section is divided into two modules, one to calculate the subsonic 
diffuser pressure recovery PRsUB another to calculate both the super- 
sonic diffuser recovery PRguP total pressure recovery to the engine 

face P^xoT* 


Subsonic Diffuser Recovery 

The empirical method of reference 2 is used. For y * 


PR 


^EF 


SUB P. 


« 1.0 - EPS 


TH 


(i.o 3 si 

I [1.0 + 0.2(M^jj)2] • J 


where 


EPS - 0. 37148 (M^) 2 - 0. 231428 (Hpjj) +0.06 
Supersonic Diffuser Recovery 

Four different options are available for calculating the supersonic 
diffuser recovery: 

(1) AIA standard ram recovery From reference 10, we have 
P. 

PR. 


■'TH 


-SUP P. 


PR 


TH 


SUP P. 


1.0 - 0.1(M - 1. 0)1-5 ; M^ > 1.0 
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(2) Military Specification 5008B — Also from reference 10, we have 
P. 


PR, 


TH 


SUP P. 


= 1.0 ; M < 1.0 


PR, 


'TH 


SUP P 


1.0 - 0.075(M - 1. 0)^-3^ ; M > 1.0 


(3) Normal shock — From reference 11, we have 
P^SUP - = 1.0 ; M < 1.0 


PR, 


SUP P 


"m / / 6 

t 'm2 + 5.0/ \7m2 - 1.0/ 


00 00 


M > 1.0 


(4) Input table of PRsUP program listing In appendix B. 

Figure 5 shows a comparison of the first three supersonic diffuser pressure 
recovery schedules described above. 

The particular total pressure recovery schedule to be used is selected by 
use of the control parameter IPR, as follows: 


IPR Code 

= 1 . 
- 2 , 
- 3, 
= 4. 


Recovery schedule 

AIA standard ram recovery — APR 
MIL Specification 5008B — APR 
normal shock — APR 
table look up 


where APR is an input incremental pressure recovery correction. If IPR is 
positive the installation effects are included. If IPR is input with a minus 
sign, the installation effects are neglected and the thrust is corrected only 
for the pressure recovery losses (l.e., IPR “ —1 gives AIA ram recovery — 

APR and no installation losses) . 


If IPR is input as a positive number, but preceded by a one (i.e., 11, 12, 
13, or 14), the installation effects are included and the subsonic diffuser 
pressure recovery is computed from the empirical results of reference 2 (see 
subsonic diffuser recovery in the previous section). Thus, IPR - 11 gives the 
AIA ram recovery multiplied by PRg^g with the installation effects included. 
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The subsonic diffuser pressure recovery is multiplied by the supersonic 
diffuser pressure recovery to give the total pressure recovery to the engine 
face. That is, 


P P 

PR « BL. X — = PR X PR 

^OT P P SUB SUP 

^TH 

Also in this module, the supersonic diffuser pressure recovery at the inlet 
design Mach number (PRdes ^DES^ multiplied by the subsonic pressure 
recovery to give the total pressure recovery to the engine face. 


INLET SIZING (MODULE SIZIN) 


This module is used to compute the inlet capture area A^. The inlet 
capture area is defined to be the total projected frontal area of the inlet, 
including the projected frontal area of the centerbody (see fig. 1). The 
inlet capture area is computed at the design Mach number, altitude, and power 
setting, and is held fixed for off-design operation. 

A useful relationship which is needed in the following development is 
the corrected airflow per unit area, which is defined to be 



WFF(M) denotes the corrected airflow per unit area (sometimes called the 
weight flow function) calculated for the Mach number specified in the 
parenthesis. For example, WFF(Mgp) means the weight flow function calculated 
for the engine face Mach number. 


Inlet Throat Area 

For external compression inlet designs with sharp lips the inlet face 
flow area is equal to the inlet throat area. The inlet throat Mach number is 
input to the program and the engine face Mach number and engine face flow area 
App are obtained from the engine description. Therefore, using conservation 
of mass between the engine face and inlet throat, the inlet throat area can be 
calculated. 


A 


TH 


|-WFF(Mef)1 

p 

*^EF 

, , 'Sp . \ent| 

[“"■(Mth)] 

P 

‘^TH 

L A A J 

L c c ^ 
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The above relation is used with the appropriate design point input values to 
calculate the design point inlet throat area. 


Inlet Capture Area 

The inlet capture area can be computed by using the conservation of mass 
relation between the inlet throat and the free-stream conditions. The inlet 
capture area is equal to the free-stream flow area (i.e., A /A « 1.0) at the 
inlet design point. Therefore, ^ 


A = A 
c o 



ADDITIVE DRAG (MODULE CDADDI) 


The additive and subsonic spillage drag computational approach follows 
Sibulkin (ref. 3). The inputs and outputs of the module are shown in 
figure 3. If the design Mach number (M^gg) is less than or equal to one, the 
bleed and bypass area ratios, as well as the additive and subsonic spill drags, 
are set equal to zero. If the design Mach number (M^^g) is greater than one, 
the following are assumed: 

1. Axis 3 mmietric cone geometry 

2. External compression inlet 

3. 20® cone half angle (THETA » 20®) can be varied internally 

4. COS X « 1.0 

5. Throat Mach number is constant at input value. 

The ratio of Aq/A^. for the engine airflow is calculated to be 


where 



p A V 

00 Q OO 

p A V 

OO Q 00 


w 

a 


p A V 
c “ 


P V 
00 00 


WFF(MJP^ 

(T 

^OO 
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The bleed and bypass area ratios are then computed from a predetermined 
schedule which can be changed if desired. The schedules are currently 


A 


c 


0,10 SFBEP 


“des V / \ 


A 


c 


SFBPP 



0.5 


where 

SFBEP = an input scale factor for the bleed flow schedule 

SFBPP = an input scale factor for the bypass flow schedule 

(Note: If the bleed and/or bypass airflow schedules are changed here, they 

must also be changed in subroutine SIZIN.) 

The ratio of A^/A^ for the inlet is computed from the engine airflow 
characteristics and the bleed, bypass, and vent airflow characteristics: 


where 


and 


Hf) 

r ' r* ' 


(1.0 + WEXWEF) 


ENG 


WEXWEF 


/^BL , \eNt\ 

W^ I A A / 

a \ c c c / 


\7ENT 

— is input (0.03 is typical) 


2 



1 

> 

n 

8 

(A 
< 


A 

c 


itt iti , , 

^ p ^ ^ 


TH 


cone , ^ o ..2 
_ - 1.0 - — 

® C 


+ c. 
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where 




A 


TH 


PR 


SUP 



W (1.0 + WEXWEFXTt...)”'^ 
WF(M^)PR3^pP^ 

00 


For < 1.0, the cone surface Mach number and cone surface pressure 
ratio are estimated, as follows: 


M » M 
cone » 



where 

pressure is 


is the cone static pressure at the throat. 


P 

cone 



The cone average 


and the cone surface pressure ratio is 


PSPIN 


P 

cone 

P 


For M^ > 1.0, the cone surface pressure coefficient can be estimated 
using an approximation presented by Lighthill (ref. 4): 
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C = cone average pressure coefficient 
S 



where 6 is the cone half angle in radians. 

The cone surface pressure ratio can be obtained from the definition of 
the pressure coefficient 


PSPIN E p /P - n X (Q/p ) + 1.0 
cone » P ^ «► 

D 


where 


Q/p - 0-7 m2 


The cone surface Mach number can be approximated by using a formulation 
of Lighthill (ref. 4) : 


M 


[j »iCp^(y + 1) + l] - MjCp (J M^Cp + 1)1 


1/2 


cone 




»Cp (Y - 1) + 1 




or, for Y * 1*4, 


M 


% + 1) 

s J 


cone 


H„[- 


(0.6m2Cp- + 1.0) - Cp_(0.35M2Cp_ + 1) 




(O./M^Cp +l)(0.1M2Cp 

s s 


+ 1 ) . 


1/2 


To complete the additive drag calculation, it is necessary to evaluate 
the subsonic spillage drag Cg. Cg is the drag of the inlet spillage that 
occurs behind a normal shock. This drag is equal to zero if the free-stream 
Mach number is subsonic. 
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Using Sibulkin’s formulation (ref* 3), we have 


C 


S 




(p/p 


cone 


- 1)P 
/ cone 


P 

00 


A 

s 


A 

c 


- cos X 


(see fig. 1) 



6 * cone half angle (see fig, 1) 


yc 


\ A 6 


) 


K - f(M ) = 0.2505M2 - 1.492625M + 2.8921 

00 OD 00 

(see ref. 3, p. 7) 

where 6 is the ratio of mass flow with supersonic flow at the inlet to the 
maximum theoretical capture area mass flow. 


Note that g is a function of ^^one^^c* ^ and, according to Sibulkin 

(ref. 3), g can be considered equivalent^ "in most cases" to the supercritical 
mass flow ratio. The supercritical mass flow ratio is presented by Barry 
(ref. 12) where g is equal to Barry’s • For the present purposes, it 

is assumed that * ° 


g = 1.0 ; for X /y < 1*2 
cone •'c 


P/P 

cone 


- »co„c'!'c - - 1-2) 

PKSPC - (7m 2 - l)/6 


for X /y > 
cone c ^ 


1,2 


For M <0.4 or A /A >1.0, 
00 o c * 


CDad - 0.0 


C3 H CDgp - 0.0 
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Figure 11 shows a comparison of additive drag coefficient as computed by 
the methods of reference 3 and by the PRINC program. Sibulkin (ref. 3) 
assumes the spike position to be a function of The PRINC method assumes 

a spike position that is a function of and throttle setting such that the 

inlet throat Mach number is a constant at the input value. 


BYPASS AND BLEED DRAGS 
(MODULE CDBYPA) 


This module computes the drag coefficients associated with the bypass 
(CDBP) and bleed (CDBL) systems. The derivation of these drag effects is the 
same; however, it is usually assumed that the pressure recovery for the bleed 
system is lower than for the bypass system. 

Two assumptions may be made for the bleed and bypass exit nozzles; 
namely, that they are either (1) sonic nozzles, with Mg = 1, or (2) fully 
expanded nozzles, with 


P^ = P 
E « 


Mg = {5 [(P^ /PJO.286 _ 


1/2 


The assumption currently used in the bleed and bypass subroutine is that 
the exit nozzles are sonic; however, if it is desired to use the fully 
expanded assumption, the changes necessary are contained in subroutine CDBYPA 
as comment cards. Also, the bleed and bypass drags consider momentum losses 
only, and do not include any drag that may be associated with the exits them- 
selves. The derivation of the governing equation for the bypass (or bleed) 
drag is discussed next. 


The thrust for the bypass (or bleed) is (see fig. 4) 


T = - P„A^) cos 

• « 

where ( )£ “ exit conditions for the bypass and m^p = m^ from continuity 
considerations . 

The thrust coefficient (based on A ) is 

c 

"bp''- ■ 'Ve * ‘'e^ ■ ®E 

S ■ -S QA~ 

c 


from reference 13, 
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A(1 + yM^) 


F _ mV + P(A) ^ 

P ■ P 

^ = X = (1 + YM2) 

where F is stream thrust , A is area, and P is static pressure* Using the 
definition of dynamic pressure, 

0 ■ I oX ■ 7 ''«X 

and using the f/p definition, the thrust coefficient can be rewritten 

P. P. 


cos e. 


(1/2) yM£ 


I yM^A 

E «> ^ » c 


E 


(1/2)A^ 




However , 


m, 


BP ^oo^P^oo “ 


P P 

I _ ^ E ^ ^ 

P P P P P^ 
00 tg t^ » E 


(11^ 


n .Ve-^Ve 




(1 + Y Mg^) 


and, from conservation of energy, 


T = T 

t. 


Using the weight flow function, which, for y “ 1*A, is 

,3 W /t” 
V a/ t 


/ 1 \ ” / J 

WFF(M) = 0.92M/ 

\1 + 0.2M^J 


Therefore, the ratio of the exit flow area to the free-stream flow area for 
the bypass (or bleed) is 


ts. (l ^ 0.2 m 2 ) "t. »- A ^ \ 

W 0.92M^ / 1 i P (it 0.2M_2 j P 

^ \1 + 0.21L2j 


and thus the thrust coefficient for the bypass (or bleed) is 
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cos 6 BP 


(1/2) yM2 c 


,11 

( P, «E VI + 0.2KJ / t-P^E P, P, P. J) 


or, rearranging terms and using the definition of gives 


D 




n.„, 

L E M_^ \l + 0 . 2 Mp. 2 / J 


r_i^«./i + o-Vvr 1 /V ^ ]\ 

1-(y/2)m2 Me VI + 0.2H.2 ) L<Pt /Pj ) ■ Vl + 0-2MjPj J J 

E 

Note: The derivation of the bleed drag coefficient is identical to the above 

derivation with the exception of the appropriate subscripts. 


It is currently assumed that 

P 

Bleed 


= 0.3 


^EF 


Bypass ^ EF 

- 0.7 p^ 


It is also assumed that both the bleed and bypass systems have sonic 
exit nozzles. 


Figures 12 and 13 show example calculations of bypass and bleed drag 
coefficients for sonic exit Mach numbers. Engine face total pressure 
recovery and bypass and bleed mass flow schedules for a study supersonic 
transport configuration from reference 14 are presented in figure 12. These 
values are used as Inputs to the PRINC module and the calculated drag coeffi- 
cients that are based on inlet capture area are shown in figure 13. The bypass 
results (fig. 13a) of reference 14, and the PRINC module calculations (dashed 
curve) are based on an exit angle of 10® and on a bypass pressure recovery 
that is assumed equal to the engine face recovery. The PRINC module results 
agree well with those of reference 14. A calculated curve from PRINC module 
that indicates the effects of bypass recovery and exit angle is also shown 
in figure 13a. PRINC module calculated bleed drag coefficients, shown in 
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figure 13b, are compared to reference 14 values for a recovery that is three- 
tenths the engine face recovery. Again, the agreement is good. Also, the 
effect of changing bleed exit angle on the PRINC module results is indicated 
in the figure. 


AUXILIARY SYSTEMS DRAG 
(MODULE CDAUXI) 


This module computes the drag coefficient (based on A ) associated with 
the auxiliary system (Cd^UX^ such as losses for cooling air ^f or various 
equipment and compartments. A description of this drag increment is given in 
reference 6. For these calculations, the total momentum is assumed lost. 


Therefore, 


"°AUX 


p a.,„v 2 
(l/2)p v2a 

00 oo n 


AUX 


where A^ux/^c ^ user input and is generally a small value on the order of 
0.005 to 0.01. 


BOUNDARY-LAYER DIVERTER DRAG 
(MODULE CDDIVI) 


This module computes the drag coefficient (based on A^,) of the nacelle/ 
airframe boundary-layer diverter system ^ diverter half angle 0^^ 

of 20° is assumed and the ratio of diverter height to boundary-layer height is 
approximately 0.5. The procedure used is to curve fit the empirical diverter 
pressure coefficients from two references: 

Reference 5, pg. 3-24, gives data atM=0.9, 1.57 and 1.97, 

Reference 6, pg. Ill, B. 4. 2, gives data at M = 2.0 and 3.0, 


The curve fit yields the following relations: 
^ 1.2 ®D \eDGE 


^®DIV H»2 20 A^ V 
. 0 . 499 ^^!^ 


for M >1.55 


for 0.95 < M < 1.55 
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- 0«8) ^ \eDGE 

(0.95 - 0.80) 20 A 

c 


0.499 ; 


for 0.80 < M <0.95 
00 


= 0.0 ; for M < 0.8 
00 

\eDGE^^c ^ input. 

Figure 14 shows a comparison of the diverter pressure coefficients com- 
puted by the PRINC module with data from the two references for various Mach 
numbers . 


BOATTAIL DRAG (MODULE CDBTA) 


The drag on the airframe back to the fuselage end point (the ’^customer 
connect” point, see fig. lb) is calculated as part of the airplane drag. 

The drag on the portion of the engine nozzle aft of this point is defined as 
the boattail drag. The boattail drag is a function of the free-stream Mach 
number, the boattail angle, and the length of the boattail. The performance 
penalty for this drag is charged to the engine performance in accordance with 
the ACSYNT bookkeeping system. The boattail drag estimation method used here 
is described in reference 2. The boattail drag coefficient is based on the 
nozzle area per engine at the "customer connect” point in reference 2; however, 
the basis is changed to the inlet capture area per engine in the program. The 
ratio of nozzle area per engine at the customer connect to inlet capture area 
per engine required for the change is 



4A 


c 


2 


The curve fit of drag coefficients based on 


A^^ (from ref. 2, fig. 41) yields 


C = 0.0102 

e 


(^) 


(1 - M^-5) 


1.4 tan g 


M 


1.53 


r / ^ \^i 


for M <0.95 

op 

for M >1.0 


For Mach numbers between 0.95 and 1.0, interpolate linearly between the 
above relations. These equations are for a nozzle pressure ratio of 2.5. 
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Values for the above equations are 

D 


ENG 






where is an input from the engine calculation. 


”cc ■ “eng 


M, 


EXIT 


NOZ 




- 1 


Izl 

2 


1/2 


(ref. 11) 


where P /P = 1/NPR which is input. 

00 t 

NOZ 




\oZ_, “ WFF(l) P 


NOZ 


W 


TH 


NOZ 


where WFF(l) is the weight flow function at M = 1.0; T , P , and W 
are input. NOZ NOZ 


\oz, 


J±L 


J[±L 


^IT 


^0Z„„ 


2(y-1) 


■ *NOZ J'^IT 
TH 


D = 
g 




^^XIT 


Assume , then 


6' = tan~^ ^ ^ in radians 


^° cc - ° 

^Sioz 

6 = 57.3 X e' in degrees 


(see ref. 11) 


To correct for nozzle pressure ratio NPR, which is an input value from the 
engine calculation, use reference 2, figure 42: 
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and 


The corrected 


ACj^ = 0 ; if NPR < 3 

6 

= 0.005 (NPR - 3) ; if NPR is between 3 and 4 

3 

= 0.01 (NPR - 4) + 0.005 ; if NPR is between 4 and 8 

3 

- 0.045 ; if NPR > 8 
3 

C is then 
°3 


C 


D 


3 



AC, 


D 


3 


To base coefficient on capture area, 



where A * ttD ^/4, as previously described. Finally, if 
to to 


< 0 , set C 


= 0 . 


BT 


BT 


Figure 15a is a plot of the PRINC module computed (based on a 

customer-connect area of 3 ft^) for a nozzle pressure ratio of 2.5 and for 
two different boattail angles. Data from reference 2 is also shoxm (symbols) 
for the same conditions, indicating the ACSYNT calculations are low for Mach 
numbers below about 0.8. 

A comparison of Boeing lightweight fighter data (ref. 2) and PRINC 
module calculations for the same nozzle (based on a reference area of 20.2 ft^) 
is shown in figure 15b. The nozzle pressure ratio for the data is not known, 
so several values are shown for the calculations. The PRINC module over- 
predicts at supersonic speeds and underpredicts at subsonic Mach numbers for 
this nozzle configuration. 


NOZZLE INTERFERENCE DRAG (MODULE ENGCDI) 


This module calculates the Interference drag on the base between multiple 
nozzle afterbodies. The procedure used is an interpolation between the curves 
(Cj) ’) of reference 2, figure 46, which have been tabularized and put into the 
program. is the interference drag coefficient between two engines for a 
nozzle pressure ratio of 2.5. The independent variables are Mach number (M^) 
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and nozzle spacing ratio S/Dg, where S is the distance between adjacent 
nozzle centerlines and Dg is the jet diameter. The value S/Dg is a user 

input to the program. The final interference drag coefficient Ct\ is bas 

INF 

on capture area per engine. For a given S/Dg and is obtained 

from the table look up for a nozzle pressure ratio of 2.5. To determine 
the final given nozzle pressure ratio and capture area, the 

following correction is applied: 


/2.5\/ ^ENG ^ \( \ 

Neng a 


X 2 X T 


where 2.5/NPR is a correction for nozzle pressure ratio and NPR is input 
to the program from the engine calculation; (NenG ” 1)/NenG ^ correction 
for number of engines since desired output is per engine and Neng input 

to the program; and Tg is gross thrust per engine for the given and 

power setting and is input to the program from the engine calculation. 


Figure 16 is a plot ACSYNT determined 


for several values of 


S/Dg compared with the data of reference 2. The graph simply shows the 
accuracy of the table look up procedures while giving an indication of the 
magnitude and variation of the results with Mach number. 


CONTROL ROUTINE (XINLET) 


This portion of the program controls the sequence of calling the various 
modules. In addition, it converts all the drag coefficients to the wing 
reference area and to the proper number of engines, since the values from the 
various modules are based on capture area per engine. 
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APPENDIX B 


MODULE LISTING 


This appendix contains the FORTRAN listing for the Propulsion Installation 
Calculation (PRINC) module for the ACSYNT program. 
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o o o o oooooo ono 


SUfc^LUTlNE LjiYPA 


SLbHjUTIN; f 0bY3A( Xrta» A^YPAC» AdLtAC»COBt»CDBP/PfiTOT#PIWF, PTD) 

C 

C Ci.MPoY!£S TnP BYPASS AND BL£cD tPFEClS 

C THE ACOniVF DRAG CALCULATION IS FOR THE TOTAL AIRFLOW 

C ff.TLRING The INLET. USING THIS BOOKKEEPING THt EFFECT OF 

C Th^ BYPASS AND BLEED MUST BF ADDED IN. 

C XE,0-ER:t stream MACH NO 

C AbYPAC*ABYPASS/AC AT FREESTREAM 

C AbLEAC*ABLELD/AC at FREESTREAM 

C CDBI. = InCREE-ENTAL drag COEF for bleed based on AC 

C CDeK=]\CKEMENTAL DRAG COEF FOR BYPASS BASED ON AC 

c rrtgt*inl:t total FkESSURF recovery to engine face 

C PINF=FRE^STPEAM static pressure (PSF) 

C PTC*FFEESTR£aM total pressure <PSF) 

C PiePPT=BYP ass total pressure recovery (.7*£NGIN£ face pres RcC) 
C PTBEPT«=BLLtD TOTAL PRESSURE RECOVERY (.3*ENGINE FACE PPcS REC) 

c 

wPT(XM)=. 7*XM+XM*(1.+.2»XM+XM)**(-3.5» 
FPT(Xh)=(i.+i.A+XM»XM)*(l.+.2*XM*XM)**(-3.5) 

PPT ( XM) = ( 1 .4.2+XM*XM )*♦ (-3. E) 

aS5U''L exit angle for bleed and bypass » 15 DEG 

>M02=XM04XM0 
CC5rt*.'V66 
PTBPPT».7+PRT0T 
PTEYa*PTBPPl+PTU 
E TBFIR.PTBYP/PINF 
IF (PTBPIN.GT.l. ) CO IQ 10 
CDBP=C. 

GO TU 20 

ASSUME THE BYPASS EXIT IS FULLY EXPANDED. FOR A SONIC 
BYPASS NJZZLl SET XMEBY (NOZZLE EXIT MACH NO) » 1. 

XMEBY = SORT (E.+ ( (PT3YP/PINF)**.2B6-i.) ) 

ASSUME 4 SCMC EXIT 

i C X Pi e b Y = 1 . 

>MtfluR= (1. + .2 + XMEBY + XMEBY) / (l. + .2*XMO*XnO) 
CDP?2»2.*(l.-((XMEBY/XM'0) + (XM£MOR**(-,5) E^CDSDE ) ) 
Cl-UP2*CD3P2+(C0S0E/(. 7+XM02 ) *xMO/ XM t BY* ( XMEMOP ♦* 3 ) ♦( 1 ./PTBPPT 
1 -<XMEMGR**(-3.S) ) ) ) 

CUBP«CUBP2*ABY?AC 
IF (CD3P.LE.0.) COBP*G. 


FOP BLEED 

2C PTBtPT».3*PFTuT 
PlBLc-PTBEPT*PTO 
pTbLPl-PTiJLf /PINF 
IF (PTBLPI.GF.l. ) GO TO 30 
CD6E*C. 

RETURN 
C 

C ASSUME THE BYPASS EXIT IS FULLY EXPANDED. FOR A SONIC 


COBYUjvI 

C03Y0C0E 

COBY0002 

CDBYOOOP 

C0BY00C5 

CDBYOCOt 

CDBY'JC 

COBYOCOe 

COBYOOOv 

COBYOCic 

COBYOCil 

CDBY0012 

CDBY0013 

CDBYOCIA 

CDBYOGIS 

CDBYOOib 

C0BYC017 

COBYOOie 

CDBY0019 

C0BY0G20 

CDBY0C21 

C0BYOO22 

CD8Y0023 

CDBYDC2A 

CDBY0025 

CDBYOC'26 

CDBY0U27 

CDBY0C28 

CDBY0029 

CDBY0C30 

C0BY0031 

COBYOC32 

C0BY0033 

C0BY0U3A 

COBYOU35 

C0BY0C36 

COBY0037 

CDDY0C38 

CDBY0-3q 

CDBY0D4G 

COBYOO<rl 

CDBY00<E2 

CDBY00A3 

C0BY0044 

CDBY0045 

C0BY0046 

C0BY0U47 

COBYOU4e 

CDBY0C49 

CDBYO^'SC 

C0BY0O51 

C0BY0052 

CDBY0053 

CDBY0U54 

COBYOL55 

CDBY0056 
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SoBkfLTlNL 


C fiLiED nJZZLF ser XM-ISZ NuZZLc tXiT iaCH) * 1. 

C xM i , + ( { PT3U-;/PU’F) **.2 86-i. ) ) 

C ASSO-^t A SCMC EXIT .'.2ZZL5 FQk Tnt BLEED 

C 

iC XMiJE^l. 

>ETM3E={l,+.3*XMEBt+XM£3E)/( 1,+. 2*XM02 ) 
CLBE>«2.*(i.-((XMErt:/KrtiT)+(xHE!iuE + *(-.5)>*COSDE)) 
CCCr2»Cl)3f:2+{C'"JSOc/( ,7*XM02)*X:’l0/XM£3£*(XrlEM0E + *3) + (l./PTaEPT 
] -(XdF-lOE**(-3.5n ) ) 

CDc’E *CD-JLi:*AoLE AC 
IF (CD’3E*LE«.-«) CDbL“j* 

IF (CDB^-.Lt.O.) CD3P»0. 

!■ E T L N 
END 


CD8YDC5 

COBYOCa 

C0BY0C5‘ 

CDBYOCc' 

CDBY0L6. 

CDBY0D6. 

CDBYODo; 

CDBYCDt- 

CDBYOCc 

CDBYOlBi 

COBYCCo' 

CDBYOCc; 

COBYOCo' 

COBYOvZv 
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o r' 


SLli^-CLTlNt P<jU3S 


V 


.-UTl N . P-.5Urti(xn,?.<Sua) PRSUC;0 

PRSUCtCC 

CL^^b^uS Tnt S'JfaSOMC OIfFUSr:k PRcSiURE a tCJ t/E=: Y PRSUOOi 

PftSUCCCJ 

X^T2 = ;^^•r*K^■T PRSUOCil/ 

I h-s* .uri^o^xMz-. 221 ‘♦^e*xMT+,t 3 prsucc:' 

PR5Lt*l.-iPS''(l.-l./U.+-.2*XMT2)**3.it PkSUOOO 

htllRN PRSUOCO 

f Lf SUOC 0 
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S Li h P I L 1 1 1 t C 0 3 T A 


C 

C 


C 

C 

c 

c 


il f<f i UTINt CuHTA(X(lG, XNPR>PTMGZ# TTNCiZ# AtNo»C 0BT» WA» AC» A£XIT> BbTA) 

X^L=f-P£c iTPLAM 1ACH NUHRcP 
XNf>=NjZZLE PRESjUPfc RATIO 
FTNn^xM,}zZL? [XIT TOTAL PPESSuk: 

TTr."jZ«\,:zZLi: EXIT TOTAL TEMPERATURE 

ALf- <01 '!l Face TLTAL APlA* so ft 

Cr.iLT = ciJATTAlL URAG PtK ENGINE, REFERENCED TU AC 

► A = £nUn£ AIki-LOr, LiJS/SEC 

AC=lNLtT CAPTURE APEA> SQ FT 

CCM'.'r’R "DRAG CURRECTIJS FACTOR FUR NPR 

Atxn = , nozzle Exit area per engine, sq ft 

V ^F ( >M J = . y2*XN / { 1 , 2»x l + XD^ + B 

XPlx = S0PT((Xi'iPR + *. 2 b 6 -i.)/,£;) 

Af.l 7T"S JkTITTfJrjZ I /PTnuZ*WA/wff( 1. » 
ANTAt«i.72e*XMEX/(l.+.2*XMcX*XHEX)**3 
At X1T= T. 

IF ( aNTAE ,GT . ) A t X I T* ANUZT / ANT AE 
AiSL RE 

A Customer connect » i.zi+aeng 
LN uZ = lia of Eng 


COBTJ. 

C06T0C 

CDBTCv 

CDBTOC 

CU8T0L 

CDBTCC 

COBTDL 

CDBTDG 

CDBTOC 

CDBTCC 

CDBTOc 

CDBTOC 

CDBTOC 

CDBTOC 

COBTOU 

COBTOl 

CDBTOC 

CDBTOO 

CDBTOC 

CDBTCC 

CDBTOC 

CDBTOC, 

CDBTCC, 


iiNii TriL MAX Exhaust diameter to customek connect diameter 

ACC » ARi.A AT CUSTOMER CONNECT POINT, SU FT 
ATKSi"1.21*AENG 

IF { AFXIT.GT. ATtSl ) AEXIT"ATESl 
DtNG"E,*SUPT<AE,NG/3.1A13')) 

DCC"1 •1'^DEnG 

DLXlT = 2.»StPT(AcXIT/3.1x,lE9) 

DEXGCC"OtXlT/DCC 
A? XIT-ATE SI 
XLFiLZ-OENG 

BITAI"ARSiN((DCC-DEXIT)/ ( 2.*XLN0Z ) ) 

IF ( JETa1.lt. C.) bETAi»0. 
btTA»S7.2'i£77'}5*B£TAI 
Tbl = l,«.*TAN(BETAI ) 

CC'.\TDI»Tbl*(l.-DEXL'CC*OEXDCCI 

IF (XM3.LE..95) C DE T» . 0 102/ ( 1, -XMO**! . 5 ) *BE T A/ 1 6 . 

IF (XMU.G'.,95) COBTl-CONTBl/XHO^^l.SB 

IF (XMO.LT..9S.OR.XMO.GT.1.) GJ TO 10 

CDBT2"CUNTBI 

CtBT*CaNTBl/.9i»**l.S3 

CDBT-CJB r + 2C.*(CDBT2-CDBT )*(XMU-.95) 

IL IF (XMC.GE.l.) C0BT-C03T1 

SET ALL CDET > (COST AT M»1.2) EUUAL TO (CDBT AT M*l,2) 

THEN CORRECT COATTAIL DRAG DUE TO XNPR VARIATIONS 

CDBTiL®C0NT6I/1.2^^1.53 
IF (CCBTiZ.LT.COBTI CDBT-CD3T12 
COkNPR*; . 

IF ( XNPR .GT .3 . .AND. XNPR.lt. A. ) CORNER- ,005 + ( XNPR-3, ) 

IF (XNPR,GE.A..AND.XNPR.U:.0.) C0RNPR-,01*(XNPR-A.)+,CC5 


CDBTCC, 

CDBTOC. 

CDBTOO. 

CDBTOO, 

CDBTOC, 

CDBTOO. 

COBTOO; 

CDBTOCi 

CDBTOC: 

CDBTOC; 

CDBTOO; 

CDBTCC- 

CDBTDC2 

C0BTCL2 

CDBT003 

CDBT-.C5 

CDBTOOa 

COBTOO<t 

CDBTOO'. 

CDBTOOA 

CDBTOC'. 

COBTOO'i 

CDBTOO'. 

COBTGCA 

COBTOC'. 

CDBTOOA 

C0BT.::5 

CDBTC05 

CD6T0C5 

CDBTOCp 

CDBT0C5 

C0BTC05 

CDBT0C5 
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SUtrKLcTIf.r COfiTA 


IF ( XNPS.GI .6. ) C0kN(>R=.0^5 CD3T0051 

Crb1*C0cT-CDP'JPR C03TCv5c 

ACC*.7b53;e*DCC*UCC COBTOObS 

CCl T=Ciin*<\CC /AC CDBT306C 

IF (CD6T.LT.:,) CDUT = 'J. CDBTOOo) 

•"'-TL'PN CD0TOfo2 

CDBTOOo3 










is 
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SLi-LUTjNl CjjIsJi 


c 

V, 

C 


SUi f.' 'UTl i.- CfiOlVI ( AWA = NG»CCO[V ) 

m DATA AT I'i = 7. IN G/D HBK PG 3.S.^.2 
ASSlf^t UlVt'eTrR HlIGHT » 

FIT jATA Al !*^.9 In 1ST AFRCl MANUAL (N/A) ■^G 3-24 

ThLlV * DIVFFTCk iNCuUOiO ANGLE 

AtiAiNo » A>-t. A Or GlViHTrP «^LOGt GIvlUcD tJY AC 

ThGI V=2^« 

CDC'Iv’* J. 

IF (aSJ. jT . .B,AN0.)(Na.Lr..45 ) CL»Jr/=>(XN0-,5) / .i5*.A«<'i*TH0IW20. 
IF (XF 3. Gt.. 95. AND.XMG.lt. CDOIV* .02495*THDTV 
IF ( <M ..G-;.I .95) CDDIV-.Ot.tTHDlV/IXfG+XKL') 

CLDl V»C0i. I ♦ AxiAENG 
F l TUi<K 
FM 


CODIUCO: 
CDOIOCu. 
CDDI DOw; 
COOIOOG' 
CDOIC' 3‘ 
CODIOCO( 
C D 0 1 0 0 D 
CODICC 5i 
CDDICOC- 
COUlGOit 
CDDIOCi . 
COD IOC IF 
CDDIOCi: 
COOIODl^ 
CDDIOCi; 
CDDIOCi.c 


32 



oor^r'. r^ooonoor) 


..I ' ■ I. ' lit' il Zl'.' 


JUT .Ni. SiZl XMEF,X*1T>P«DCS#AC>Xr!L;fi/PhSiJB#SFBEP» 

C i;*. 'rT i N. TU SIZt I TLETS 

C 

V ; ; 1^^.' iL7tL 6K‘, bYPASS SCTEJULEi FGH iNlET 

^ >•' TL Tn^T TH:SE SCHEDULES MUST PT COMPATIBLE kITH THE 

^ ^LLt;, .i Hi PYPAS:; SCHuJULES IN SU3-lJUTll.t CDAoDi 

C Am * FACE FLCiW A^EA, FT + FT 

XF;F » _ '<oifG FACl FLOrt MACH Xhj 

XM « THFIAI MACh NT 
PK « SUP^fCNlC jIFF. P,R, 

P^SI J « juFS'NIC ulFF. P.P, 

ATU = D.iU^ THSJaT FLUV APcA» FT*FT 
.T * THROAT ^klA, FT^FT 

At » JNL'm LfMTUPt AFEA, FT^'FT 

= i\l:t Design mach 

fpT' T » TIT AL P-cSS -(.C. T j tF 
Zi-bFP = SCALE cACTGk p CJ < INuET BLEED D^aG 
St-^P^ « SCALr FaCTD< FUP IN'LiT oYPASS DPAG 

fr-lF^i*.i*(XhDES/2.)**’*SFBE'> 

AvLaCL*.j 3 
tf F Xr'T = WFF ( XMT) 

JTL=a>.F*( WFF (XMEF) / wFFXMT)*PPSJ6 + (l.+AVtACD) 
At»ATU*(*trFXflT/wFF(XMjtS))*P90£S*(i.+AijLACu) 

If (fC.L-.ATO) AC»ATU 
F t TLFN 
L FD 


SIZID'-'J 
SIZIO-J 
SI zi jcc 

sizio:.- 

SIZIOOv. 
SI zio: D 
SI ZIOUv 
SIZIOOO 
SIZI jC v' 
SIZIOOl! 
SI ZIOCa 
SI ZIOC*. 
SlZICC'l. 
SIZICl i 

sizioui; 

SI ZI»ui! 
SIZI u* 
s I z I c : . < 
SIZlOul- 
si ZIOUJC 
SIZIuC Z. 
SIZIOO^; 

siziotz; 
SI ZI002F 
SIZIuCZI 
SIZI J02t 
SIZIC02T 
SIZIOC2E 
SIZIOC29 
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C 

C 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


1 < 


2 C 


Jl 




5ct'K(jtlXNL fNbCDI(XMU^nN^SGDG^tDI>FlP»0#JiCf^4'iP^) ENGCOvW 

LIr’FNSiGN ID55 (^)>CGp 5(^), SD75< ^}»Cu75 ( 4 ) , SD 8 5 ( 4 ) , C C 0 5 ( 4) > SD95 ( ^ EN GCOO-^ 


2 CL9;(<.),3C 1w(^)#CD:C(4) 
L/'T* 

I>031//1«> 1 
SJW/I •# ? 


Data 
r A T A 
L A 1 A 
b/TA 

c r? ] * ‘V • 

H ( Xm'1,LT.1,2) go ro 
IF (Xf'^.j«Lr»l*b) GL* Tu 


CD53/.1}5,, j: 2^/ 

. 3 > 1 . 3 . 5 / , C 0 75 / • G 35, • J i 3^ C b / 

. 7f 2.1, 3, 5/>CD35/ ,036# • ^Cb/ 

• 2i 2.6, 3. 5/, CD95 / .UG7, .')37:>,2^ .025/ 

• 3f2.7,3.3/,CDli. / .3J0,.‘^6?,«*3b,.''28/ 


. . 0 . ) 
3 J 
ID 


RETjFn 


XFC GKCAirP. 1HAN Uk tIUUAL TO 1.8 


ENGC JOG 
EN GC GG ^ ‘ 
ENGCOlv 
GNGCOOG 
tN GC DC C* 
ENbC'Jvt.^ 
EN GC u G<^ ' 
ENGC OC X*. 
tNGC DOl ; 
ENGCOOx; 
ENGCOCi: 
EN GC OG 1 ‘ 


CL 12 )C^-j 

Gl TG IGO 


ENGCOCI: 

ENGCCGl^ 

engcolx'; 


X N L* 3 L r t Fi i • 2 A N L 1.3 

Ih (>FJ.bT.l,5) Gu TO 2> 
f i'*(X^J--1.2)/.3 
cm: = . 
ct Tj 

Lr=( xr^'2-1 .5) /.3 
CDU = .DC7;j-.u025^^D!^ 

CG Ti. 10 j 

>i"r LESS ThAN 1.2, TABL- LiJuKUP KEOolKtb 
DM Tiu* D • 

r:rTi>V'L»D. 

Jf (XA-J.GT,!,) GU TU eO 
IF <x:‘iL).Gr..95) GO TU 7D 
JF ( XF*a.GT. .65) GO TO 6> 

If (X^U.GT..75) go T3 bO 
IH ( ^MJ.GT. .:>5) GO TJ 43 

X^p[; BcTXEEN 0. AND .55 

CALL taint ( 505^ f CD55, SGDG, C OU, ^ , 1 , 1'ilRR , ) 

C b L « ^ • 

DM AH* *55 

G( Tl, 90 

XF'D 8ETwLEN .55 AND .75 

CALL taint (SC75,CD75,SGDG,CDU,4,x»NEKR,DM 0NU) 
CALI T AINT ( SL‘55f C0 55, SOOG , C DL , 4, 1, N E RR > DMONL ) 
[MTA3».2 
CF *XNC-.55 
CD TO 90 

.75 AND .35 


ENGCOClc 

tNGC 

ENGC002:. 
ENGC002I 
E N GC C U2 2 
ENGCCC2; 
ENGCOOZ^ 
ENGCOGZf 
ENGC002t 
ENGCCt27 
ENGCOOZfc 
ENGC002^ 
ENGC003C 
ENGC0031 
ENGCOC32 
ENGCOD32 
ENGC0034 
ENGCGO 35 
ENGC0O36 
ENGC0037 
ENGCC03b 
ENGC003<^ 
ENGC0O4i- 
ENGCOGhI 
ENGC0042 
ENGC0C43 
ENGC0G44 
ENGC0O45 
ENGC0046 
eNGC0047 
ENGC0O48 
ENGC0C49 
ENGCDG5C 
ENGC0051 
ENGC0052 
ENGCCIC52 
ENGC 00 5^ 
ENGC0055 
ENGC 0056 
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StttKLLlIKt tr<GCLl 


3' CALL TAINT ( SL35>C0bi» SODG>CDU»A#l»N£kS»l?fHjNU) 
CALL TAlNT(St75»CO7i;,SQ0G»CDL#A»l»NERi«#DHONL ) 
DHAB«,: 

LT' = XNC-.73 
to ID -).- 

Xf'L- BiT«»FtN «85 AND «95 

ti Call TA1nT( S095,ct)^3* SOOG»CDU#<.»i»NERR»DMiJNu) 
CALL taint (SD85»C085,SO0GjC 0L»A» i»N£RP»DMJNL ) 
CT TA3». 1 
PM*XM0-.bi 
tl Tu 9C 

Xl^r TtTWtET .95 AND *.0 

7C CALL TAlNT(S01O,CL10.SOL)G»CDU»A,i,NERR»OKUNU> 
CALL Ta1nT(SD95»CD95#S0DG,C0L»^» 1»NERR.DM0NU 
DN7AB» . 03 
CN»Xf’G -.95 
Gu TO 9C 

Xi-'li otTwcEK i.O AMO 1.2 

BL CDLi».C19 

CALL T AlNT ( S01O,CD10»SO0G»C0L» A# l»NERR»DMaNL ) 
OrTAB - .2 
C(' = XMD“1. 

vC CDIl»CDL+(CCU-CDL)»Ort/OiTAB 

DElERMNt CDI FOR T^£ ENGINES 

IOC C0XT=(2.3/XNPR )*CDI1*2,*FIP/(0»AC» 

CLI*(EN-1. )*CDIT/EN 

RETURN 

FNL> 


ENGC 0 C 57 
ENGC0C58 
ENGC Jt' ;<9 
ENGC0060 
ENGC0C61 
ENGC0C62 
ENGC0C&3 
ENGC0C6A 
ENGC0065 
EN GC OG 6 & 
EnGC b 7 
ENGCGC66 
ENGC0069 
ENGC0G7C 
ENGCCC71 
ENGC0072 
ENGCUG73 
ENGCD07A 
ENGC0075 
ENGC9C76 
ENGCGC 77 
ENGC0C76 
EXGC0G79 
ENGC0060 
ENGCjOB 1 
ENGC0G62 
ENGCOC63 
EN6COOb<I 
ENGC0C85 
ENGC0086 
EN GC 00 97 
ENGCOG80 
ENGC0039 
ENGCG090 
ENGC0091 
ENGC0092 
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iLBt-Ouri M: CDADDI { AT> AC» X»',J,pP, XMr» XMO t S # P TU» T TO » CD A D» C DA 0 S » 

I A Y r A C # A ■ 5 L t A C » A :j A C » W A » S F 3 e P # S F B P ;» ) 

C 

C AT » IWl^T ThRuAT FLTw AREA* FT+FT 

C AC * iNL-r CAPTURE ARtA» FT*FT 

C X»r = FREE STFcAh ^aCH 

C RF » SUR-^SONIC DIFFUScR TOTAL PRESiU^t PfcCUVcRY 

C XFT * Ii'tLtT ThRJAT HaCH NO 

C X'U!:'^ « INLET DESIcu MACH NO 

C tli; = FREE STREAN total PRtSSURt 

C TTl = Fhlt STREAM TOTAL TEMP 

C CDaD » SOPiFSONIC SPILL AOUITIVL DRAG BASED ON AC 

C CLADS - SUBSONIC SPILL ADDITIVE DRAG 3ASFO E).\ AC 

C AC'AC • A'J/AC 

C AFLLaC = AELt/AC 

C AETPaC = AEYR/AC 

C LA * ENGINt AlRFLOkf LBS/SEC 

C Xr'CC’NE » CONt SURFACc MACH NUMBER 

C CFCS ■ Cjv(L SURFACE PRESSURE RECOVERY 

C PSPIN . PSTATIC ON CCJNE/PSTATIC FREE GTkEAM 

PNSPC * static press, ratio across N.S. at C.-NE SURF MACH 

FRTC1(XM)»n, + ,2tXM*XM)+*(-3.5) 

V.FF(XM)s,^2*XM/(1, + ,2*XM + XMI**3*,3D:1 
c 1 NF I N ( X M ) = , D SO 5; + ,7*XM*XM 
FF.SPCni X‘il) . ( 7. ♦XMl^XMl-l, ) / 6. 

LT FXr,0* .<r F ( XMO) 

ACAC G*WA*SQF T(TTa) / (»|FFXMQ + PTC*AC) 

ALiAC * AOaCG 
CS»0. 

xfu: » XM j* xmo 

IF ( XMDES.gT.I. » GC TO 10 
CD/iL*C. 

CDADS=J. 

ADYPaC*^. 
aelfac*'. 

R i I U' Fi 
IL GAR*i,A 
ThhT A*2C. 

TLUTA1»ThET A/57.29S77V5 
CcSl AM»1 . 

tkl next Cards are the bleed amo bypass schldules for the inlet 
these Should de made compatible with the inlct d'^sign point 

RtLU;S IN SL'SGUTINE SIZIN 

ADLF. AC«.l+SF6EPR(XMDES/3.) * + 3*( XKO-1.) /( XMDtS-1. » 
AbYP»C«.D*SFBPP*( I.-AOACG) 

IF ( XRJ.GT. 1. ) GO TO 20 
CUAC«C . 

AbLEAC*C 

XNCCM »{XMO + XMT) /2. 

PSP]H»1. 

^0 IF ( AOACG. GT. ,97> AaYPAC*0. 

AViAC* .C3 

WCaP=AC*WFFXMO*PTO/SORT(TTO) 


COADOOU j 

CDaOJCCE 

CDADOCUE 

CDADOCGR 

CDAD JCOf 

CDAOScOt 

COAD0007 

COaDOOJc 

CDADOCOS 

CDaDOC j, ^ 

CDAOC’Gll 

CDAD0C12 

CDaOOcIj 

COAOUCIA 

CDAOOCl; 

CDAD0016 

CDADOwi? 

CD ADCf.ie 

C0AD2C19 

C0A00C2O 

CDA00C21 

CD ADO- 2 2 

COADOC23 

CDa0002h 

CDAD0l2J) 

CDADUC26 

COADOC27 

CDAD0028 

CDAOOU29 

CDA0-D3C 

C0AD0031 

CDAD0032 

CDADOC33 

CDAD003<» 

C0ADGO3S 

CDAD0036 

CDADC037 

CDA00C36 

CDADG039 

CDAOOCGO 

CDADOOMl 

C0AD0C^2 

CDAOOO<*3 

CDAD00^4 

CD AOOG4 5 

C0A00C46 

CDA00CXF7 

CDA00C48 

CDAD0DM9 

CDAD005C 

CDADOOSl 

CDAOOC52 

CDADOOi»3 

CDADOCSR 

CDA00C55 

CDAD0C56 
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ri.'HKfLUf-L iJMjOl 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

t 


c 

c 

c 

c 

c 


V.', Xir EF a ( A n 1 AC -t-Ad YR \ C + AVE AC ) *wC Ar / XA 

CJ AD0C3 

Al AC = i 0AC3* ( : .4 -^.-x,EF ) 

C 0 AO DC :> 1 


CDaD0C3‘ 

If L i T G: UAL TRY 

COAOOCoi 


COA0OO6: 

T f'i T V T a A F F ( X 'Vn ♦ ‘5R ♦ P T G 

C D AD wyv Oy 

/.T=*A+(l. + v,tx«tE)*SQPT{TTO)/RHJTVT 

coADOCo: 

Ii (VYiu^,AC) ATa,^9<i<AC 

Cu ADOCo^ 

<*S = «L- AT*Cf..SL AM 

C 0 A 0 0 w & ! 

t't[n=3wRT(A3/3.1‘ilci) 

CO ADvG Cc 

y/LLNEak.: jNr/(TAN{ thc iad ) 

COAOCOc* 

YC af .R T ( AC / 3 . *.A1E ) 

CO AOOODt 

XLt YC=XAC uNF /YC 

COAOOub^ 


C0400C7: 


NlTL r.-in The INLlT GEOMETtiY HAS ^EtN SPECIFIED PY ASSUMING THE 
THF'CiAT MACh number AND THE ANGLc OF THF GONE Ok 54MP 

AT = iNLcT Throat now area# ft*ft 

AS = frontal area CF inlet C/B at inlet T.TRuAT/ ft*ft 

wClNc - RACiUS .if CuNr AT ^^ILFT THRCAT# FT 

>ACDNt » OISTtNCE tHuM CONE TIP TO INLET THROAT, FT 

Yc « kA.jiJS cOFRESR J nDING 10 INLET LAPTUkL AREA, FT 

>CUYC a XCCNE/YC 

K’R RAtP OF CONE static PPcSSUkE USfc AVERAGE OF 
FkLL stream and throat static PKcSSuRE 

hiR . =Pk vkkTGK xMTI/PPTOTtXMO) 

CPC3 = .5R( - lK-1. ) / CINPIN(XMJ) 

PSP JN» ( t .p +.,)/?, 

IT ( X.MO.L J.l. ) 30 TO 30 

Ct.NE SURFACE PRESSURE RATIO 


C0A00C7J 

C0A00O7E 

CDA0uu7; 

CDA00C7‘ 

C0ADLC7t 

C0A0CC7c 

COADUL77 

C0ADJG7C 

cdadg:7s 

CDAOOCdC 

CDADOOol 

CDAOOCe^ 

C0AD0083 

CDAD'JUdA 

CDADCOdS 

CDAD0C66 

CDA00&87 

C0A00088 


AMLil = X.M01-l. 

S'->MC1 a s:jt -| ( XMOli ) 

AH/=ALOo( E./( SaXMOl+THETAI n 

THtTAAaTHrT/i**-A 

A.=THEr41++?*(2.+ALG-l.) 

A?*3,*X.’'.l111 + ThETa^*alG*AlG 

A3 = -(S.*X.Xri-l . )+TH£TAA*ALG 

Attain' TA<,a( 1 3.*XM01/‘i,.+ .5 + 2.A*XMLH*XHQi/X.’i011 ) 

CFCS=Al+i. 3 RA 3 + AA 
IF (CPCS. 3 T..R) CPUS *. 9 
PSPINaCPLS’ROINPlNlXilCD + i, 

3C It (XXJ.GT..G> GO T'j AO 
LL AD = '; . 

CL’ADSa 0. 

RLTUkN 

Av CDAa3,/(GAC ♦XhOl) 

CL AL = AT/AC*FiPC"»(GAM + XMT*XMT*i. )*CaSLAM 
CoAOa(AS/4C)*PSPIN-l,-AOAC*GAM+Xi10l 

Tht aPoRjX for ThF cone SURFACE PRcSSUPt CJEFF INTRODUCES AN 
ft-Pl.R Ik the additive DRAG CAL WHICH CAUSES THE ADDITIVE DRAG 
Tl be non zero at AUAC»i.J THE NEXT CARDS iNTRCPUCt A 
CUTR^CTlJN TO THE ADDITIVE DRAG TO COMPENSATE FOR THIS ERRDP 


CDAD0U09 
C0ADOC90 
C0A00O91 
CDAD0C93 
COADC 093 
CDADOU7A 
CDADOC95 
CDAD0C^6 
CDA0at97 
CDAD0C98 
C0AD0O99 
CDADOlOO 
COAOGIOl 
CDAD0102 
CUAD01D3 
CDAOOlOA 
COADO105 
C0A0 d1j 6 
CDADOJ.C7 
CDAD0108 
CDAD0109 
COADOllO 
CDADOlli 
CDAD0112 
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SUb^LL'l 1\1 C->\D:iI 


c 

All* AC 

AUL = i-,rf ( Xf'T ) /PS+ADl 

/ 31* TDl, « CuSL 4M 
C[;tr;:i;i,*(A3_/AC)*^SPlN-l,-GAM*Xi'101 
CGanil»4T0r /AC + PIPC* (GA1*X'1T*XM + 1. J *CuSL« -I 
Ci.CGr = -CG4* (L JAD]l + Cu4D21 ) 

U (XMrJ.LI.l.) G'j to ST 
XK = .21 GS♦X^ fJl-l , A92 £j 2S*XH0 + 2 .o 921 
C 

f'-TA * Aj/aC FJP MiaMUM S JP£kC1<ITICaL SPILlAGl 
I' tlA = FCN(XCnrC»Xf'0> THETA) 

I -- Ph-'i-'ir uS3 APPAIO* ^ALUC FD« 20 ObG COHc AND XMD«1.A 


F t T A ® i . 

If (xCuYC.Cr ...?> StTA*l.-(XC0YC-i,«:)/1.55 
U (.3tT4.Lt....) 3FTA*. TfOCl 
TAfiTH=TAr)( THtTAl ) 

XIYC = XX*(1.-AlAC/Bl-TA) 

AY=At*(SJ.^T(AS/AC)-XLYC*TANTH)**2 

CL'AE SukFACL MACH NU*1BE.P 

3] sX'A: l + CPC S 

t;t = .c>-)l + l . 

tr=. . 

E A = . 7*31 + ; . 

XFlCNt-X.-lO + SOETI (b2-CPCS + B3) / (3<t*iiS) ) 

IF ( Xf'.CCN- .LT.l. ) X'lCONi:*!. 

FfJSPC=?KjPCC (XrCQNL) 

C3.= (2./(GAf'*XH01))*(( AS-AY ) / AC » ♦( Pf(SPC-l. ) ♦PSPIN 
:>L Ct At*CUA*(CCAi + COA?)+CDCOP 

VLNG AND ANDERSON FORMULATION OF ADDITIVE DRAG 
PDF NASA TN 0-7AA5 

Vlvf *( .2*XNT/XM0)*SJftT(n.+.2*XM01)/(l.+XMT*XKT)) 

(-PST..* (PlPL-1, )/OlNPlN(XMO) 

CDaDZ = 2.*4CiAC + ( V1VO*CO$LAM-1. )+CPST1*AT*CU5LAH/AC + CPCS*AS / AC 

ccA(;z«cr>A.)z+CDCtjR 

clmI’S»cs 

IF (ADAC.LT.l.) 30 TO oC 
CLAD«C. 

CCADS*T. 

Pc TUX N 

ot cF (LDAD.LT.r. ) COAO^O. 

IF (CDADS.LT.C. ) CDA0S«3. 

F tTLRN 
LND 


C0AD0113 
COADDj.1 A 

CDAO.niS 

COAOOllfc 

COADOii? 

CDAOG^i e 

CDAD0U9 

CDA00120 

CDADC12i 

CJA00122 

CDADC123 

C0ADO12X 

C0ADO125 

C0AD0I26 

CDA00127 

COADOiEB 

CDAD012 

CDAD0130 

coaog:31 

C0AD0132 

CDAD0133 

CDAD013A 

CDAD013S 

C0A0i)136 

CUAD0137 

CD AD013B 

CDA0013<y 

CDADOxAO 

CDAOOlAl 

CDAD01A2 

CDAD01A3 

CDA031AA 

CDADOIhS 

CDA001A6 

CDADU1A7 

CDAD0148 

CDADDit9 

CDA00150 

CDADOIS I 

C0AD01S2 

CDAD01d3 

COADOl jA 

CDAD015S 

CDAD0156 

CDAD01S7 

CDADJ158 

CDAD01S9 

CDAD0160 

CDAD0A61 

CDAD0162 
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tSxCtllNE 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


XP Li: r( Pru^TTn,?TNJ2f TFNOZ# wACf PImF#A?nG# FRSUb# 

1 (■ K I T 0 T # A T / K L Y Z ) 

CLh' .N/KKCC'^ / Z1 ( 3) >tM,22( i6)»S«ING»Z3(2)»CDlNSP»Z4(66) 
C[lN^N/'>f^1C/XHEF,Zt)(15),XK3*Z6(7).AtFW7(9) 

Ci f ;N/Pv Jif i>/xiPR,D;LP^» IPk»XiT» XHDES»AUA :fiG, A w AENG / SC3G / 

: (6) , XPRI (5), y1,C0aFTP,Y 2,C£)INLP# Y3,CDIP»CDBTP,CDDH/P,CDAUXP, XINLviC'Cfc 

2 ClhPf.C!JJLP,CDAjS?>CjA0P,A0AC»YA(3)»PCtFAC#SHNSP»SFADP»SFADSP> XlNLu:-C7 


XiNLOJCi 
XINL<jO'^Z 
XINLOL03 
XINL0004 
XiNLVi.C P 


■ Sf BcP,SF3PP,SF AjXP.SF0rVP,SFIP#SF3TP» A&LFaC»A3YPaC»BFTA 
C i r f JN / 3 T J ^ A G / Y 5 ( A ) » AC » A N u Z » Yo ( 1 H ) 


FlcTIKF 
Xi'L 
L ■ L P K 

IFF’ 

AFF 
Xt LF 
XU 
> FO 1 S 
PT'i 
TU 

AUAc N6 

AVa- NG 

XKFk 

FIfCZ 

TU,uZ 

5CLG 

b 

FJF 

t> AC 

UNF 

ALf>iG 

FFSU-3 

Fk 

Pf TCI 

AT 

AC 

CtALP 

CLAC SP 

CLbtP 

CC'BP'^ 

CLAUXP 

CCDI VP 

CCFT? 

CDIf 

CL iKLP 

CCAFTP 

CLINSP 

ADAC 

ABYPAC 

ABLL AC 

KEYZ 

fcdfac 

Sf<<<< 

ANOZ 


TO CaLCULaIL install AT 1D.N LOSSES 
» FPEE SIRcAii MACH NUMBER 

» INCRtFiNTAL PR-SSURfc RECOVt»Y kt.JUCTION IF.POT 
AS A POSITIVE NUMBER 
PPLSSJRL RtCOVEFY CUGE 
EFGINt FACE FLOW ARtA# FT+FT 
1-FGlNc Face MACh NUMBER 

ihkdat mach number 

INLET DESIGN MACH NUMBER 
FF.EE STRlAM total PRESSURE 
FFtL STREAM TITAL TEMPERATURE 
AUXILIARY AREA OVER AC 
LlVcRT;P WiUGE AREA OVER AC 
NOZZLE PPcSSUPE PATIO 
NOZZLE EXIT TOTAL PRcSSURu 
NOZZLE EXIT TOTAL TEMPERATURE (k) 

NOZZLE SPACING OVER JfT OIAMtlER 

FREE STRLAM DYNAMIC PRESSURE 

ENGINc GROSS THRUST PtR ENGINE 

ENGINE CORRECTED WEIGHT FLOw 

FREE STREAM STATIC PRESSURE 

ENGINE FACE TOTAL APtA FT*FT 

INLET SUcSOMC DIFFUStR PRESSURE ktCOVEPY 

INLET SUPERSONIC DIFFUSER PRESSURE RECOVERY 

INLET TOTAL PRESSURE RcCOVcRY TO ENGINE (-PF+PrSUB) 

INLET THROAT FLOW AREA (FT*FT) 

IRLtT CAFTURE area (FT+FT) 

INLtT ADDITIVE DRAG PER A/C BASED ON WlNG AREA 
inlet subsonic SPILL DRAG (PER a/C) SWING REF 
INLET 3LLcD DRAG (PER A/C) SwlNG REF 
INLET 3YPASS DRAG (PER A/C) SwING REF 
INLET AUXILARY AIR DRAG (PER A/C) SwiNG REF 
INLET OIVtRTcR DRAG (PER a/C) iwlNG REF 
NOZZLE BOATTAIL DRAG (PER A/C) SWING RFF 


Base drag for space between engines 
inlet total drag per a/C Swing ref 

AFT END TOTAL DRAG PER A/C SwiNG REF 
TOTAL PROP INSTALLIATIUN DRAG PER A/C 
FFEESTFcAM FLOW AFSA/AC 
FFFESTPEAM flow area for BYPASS/AC 
FRtcSTREAM FLOW AREA FOR INLtl ELEEO/ 


(PEP A/C) SWING 


SWING REF 


1» INLcT DRAG-PR COMPUTED; -2x 
SCALE FACTOR FOR INLET CAPTURE 
SCALE FACTuk FOR VARIOUS DRAGS 
<<<< CORRESPONDS TO INSP» ADP# 
NCZZLE EXIT AREA> FT*FT 


AC 

AFT END EFFECTS 
AREA 

ETC. 


COMPUTED 


XINLOOOe 
XiNL Ot ^9 
XINL jC'ID 
XINLOCII 
XINL3L12 
XINLwUli 
XINLOUl^ 
XiNLOCiS 
XINLUC16 
XINLOOi7 
XiNLOClb 
XINLU0I9 
XINL0C2C 
XINLOuEl 
XINLvD22 
XINL0C23 
XINLUC2R 
XINLUC23 
XINL0U26 
XINL0027 
XINLOOZe 
XINLCC29 
XINL0C30 
XiNLOOi 1 
XINL0032 
XINLD033 
XINL003A 
XINLO035 
XINL0G36 
XINLOL37 
XINLOC38 
XINL0039 
XINLCDRC 
XINLDO^l 
XINL00^2 
XINL00M3 
XINLCC<i^ 
REFXINLOC45 
XINL0C46 
XINL0047 
XINLOO<»8 
XINLOt^9 
XINL0C50 
XiNLOOSi 
XINL0052 
XINLOD53 
xiNL00;j4 
XINLOU55 
XINLOD56 
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n o o o 


1 iiit X 1 NL rT 


C 


c 

c 

c 

c 


C L* A L j P “ .'' • 
c uu p r : , 

C L- P P - 3 « 

Cl / u X P = . . 
c u - 1 V p * ,? . 
c li T p = . 

Ci-iP = C. 

CuiNLP* .. 

CLAP fp* : • 

CClNiP=>, 

L-.7»xv,-:.«x.^r*PiNp 

(\EYZ.:-c..M GO TO ;o 
r P 5 1 - i • 

Tra INLtI IS ASSUMED TO HA\/c G?: U Tri Y THAT CAN BF VARIfcO IN 
SLCH t MAS\Ek TC KtEP THE THROAT HACH NuMdSR AT THE INPUT 
VALU:-. THE INLET MCUcL IS (FOr! THE PRES-NT) MIR AN LXTcRNAL 
CLhPKESSION aXISYMETRIC INLET. IF IT IS DESIRED TU NOOIFY THE C0DEXINL0L7 


TC hA^:DL^ A FIXED INLlT THE TH.HUAT FLJw AREA SriGULD BE SET XINLCC7 

FwUAL TO THE DESIGN VALJE AND THE SUBSGNIC DIFFUSER PRESSURE XINLjU7 

rtCLiYERY SHOULD BE COMPUTED FROM CON F i NU IT Y « XINL0C7 

XINLOOB 

If (IPK.oT.iO) CALL PRS'JBS(XNT,PRSUB) XINLOOB 

CALL PRlNUXMG>DtLPR»PR»IPR»P«rUT»?RSUBfXM?RI»XPRI»PRDES#XMDES) XiNLOOB 
Ir (IPR.LT.C) RcTUkN X1NLOC3 

RA-.Ui07:.D*hAC*PTD + PRT0T/SQRT(TTU) XINL0C8 

Call S I a 1N< a- F, Xrtf F, XMT» PRD£S» AC» XMOESf PRSUiJ» SF6EP > XINLOOB 

CALL CDADDI (AT. AC» XMO>PR>XMT, XMOFS.pt J#TTJ» CO AO, CDAOS.ABYP AC » XINLDOo 

I iPLEACAOAC.WA.SFuEP.SFbPP) XINLOCB 

CALL CDBYPA( XHG, A0 YP AC , ABLE AC ,C C 3E , C DBP, PkTGT, P INF , PTO ) XINLOOB 

CALL CDAUXK AUAENG.COAUX) XlNLOUo 

CALL C Julvl ( XMO, AWAENU.CDuIV) XINLOOS* 

IF (KtYZ.tC.l) RETURN XINL0C9 

It CALL COBTA(XhG,XNPR,PTNJZ,TTNGZ,AENG,CDbT,WA,AC,ANOZ,BETA) XINL009 

Call ENGCDI ( XMO, EN,S00G,C0I,FIP,Q,AC,XNPR) XINL009 

FATlD*4C«ENAPCDFAC/S*lNr, XINLOOB 

CLADP*CDAD*RATI C+SF AD? XINLCC^ 

CLACSP=«CDACS*PATiri*SFAOSP XINL009 

CDtjEP = CDJE+PATIO XINL009 

CCBPP-COBAAF ATIG XINL009 

CLAUXP-C JALX*RATIO*SFAUXP XINLQ09 

CDDIVP»COi)IV*RATIQ*SFOIwP XINLOl J 

tui;TP = C J8T*F aTIOASFBTP XINLOlO 

CD1P*CDI*R ATIj+SFIP XINLOlO 

CLiNLF»CDADP+CDAOSP+CDBEP+COBPP+CGAUXP+CDOIVP XINLOlO 

CCAf TP-CDBTP+CDI? XINLOID 

OdINSP*(CD1MP-fCDaFTP)*SFINSP XINLOlO 

return XINLOlO 

END XINLOi-o 


XINL0C5 

XlNLUi^t 

xinlocd 

XINLGOb 

X I NL ■jL o 
XINLOC c 
XiNLOCo 
XINL0C6 
XI NL jC c 
XiNLOOc 
X I NL OC 6 
XINLOC'6 
XINLCOo 
XINL0C7 
XINL 0 C 7 
XINL007 
XINLOe? 
XINL0C7 

XI NL De7 
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il'Sl- f CTl 


C ) A Ll X 1 


.’L'TlMz CLAUXK AUAtNG^COAUX) CDAUOOu 

t CDALCU'J 

C Ml-: ^.'lT;:^^AL AERDJYNAMICS MANLAl» NAR, K-? 7-£A SE-C 7.8 CDaUw-CU 

C CI;Lx - AuXlLlAr!Y SYSTEMS DRAG COAUCC 

(. ALAi-;C =■ AREA OE AUXILIARY SYSTEMS OIYILED BY AC CUAUCC 

C CDAUC-Cv 

LLAUA*2.*\UAEMG CDAJOCU 

FlUkN CDA'JOCO 

tfu coAUj:-'' 




41 


C) c. 



o o o o r: 


51ckCL 7IK P^l-iL 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


Lt:.F Pt^INL(Xhu,DcLPC^ tPk>iPk >PnrL}T,P. SU^#XMPi>I#XPRl,P>^L-S# 

J L/ 1 S ) 

bUt-K .UTI... Ti.' CC:«)PUTn THt INLET TJTAL Pt'^rbSUKE -JECOVEPY 
l^)■ = P"tiSLii£ PECOi/LkY 3PANCH CJOE 
» 1» UA STANOaPD 
* Z • f' i L SPEC :■ 'j : 4 B 
= 3» n:.F''4L shtck 
= TAtLlr LDUK-'JP# PP VS .1ACH 
x^>, * P-icr SFrdA.'l KACH '(UN3EP 

LlLrP = I'.CPi-'IcNTAL PPrSSbPS KdCOVdPY P^OUCTIu!'. 

l iHT AS A PCSITIVl NUMfiSK 
PMl,I =■ TjTAL Pk Tb inlet FACE 

L JAtiSI J J AFPP1(6),XPPI(6) 

I f t DUh* : 

X f !. 2 = > i : ♦ X L, 

X A’ L 2 = X /, C c 5 + > Ai C t S 
if ( I F R . J . . r ) b 0 TO 10 
1 Pf OoA' * I PR 
If f= = 1 A3i ( I PR ) 

* ' i F ^ 1 P P • 4. • 1 ' ) b J T 3 C ~ 

I* H CUM » I op 
JFk*1PR-.} 

PC Ct Tv { 3 j»AC, 50» 7C)f IPR 

Alt STanL'.'.hP 

3C rK=j. 

P r tE S = 1 . 

lA (xfo.^r .1 . ) PK»i .-,!♦( xMO-i. )**i, 5 
IF ( ffA!D;v,cT.i. ) PFOES*!.-. :*(XM0ES-1. 

G L T J 6 - 

ML SPEC i>CCit 


PP iNvv vl 
Pk INyOiZ 
PRIflOJCS 
PR IN ,'v> H 
PRINOC OL 
PP INObJo 
PR INCCt 7 
PKINS'.',8 
PRIN w G 0 •» 
PRINOOiO 
PR IN Or *1 
PR Ink 1 2 
PR INCCl 3 
PR INOOlt, 
PRINCvli) 
PR IN0016 
PRINCJ 17 
PRINfCie 
PR IN0C19 
PR INGt 2C 
PR INOO^i 
PRlN0t22 
PR IN0C2 3 
PRINJC2R 
PRINCU21* 
PRIN0C26 
PRINOC27 
PR IN0G2R 
PRIN0029 
PRINCf '^30 
PR1N0L31 
PR IN0032 
PR INOt 3 3 
PRIN003A 
PR1NC4.35 


FF= 1 . 

Ffr;ES»l. 

IF (XPj.oT.l.) PR»1.-.C7S*( XrtO-i. 35 
IF (XMJEi.CT.l.) PR0ES=L.-.C7i+(X*l0ES-l.>**1.3S 

GO 10 fc', 

NLRAAL SHJCK 


it Pk = 2 . 

PROtS'l. 

IF (XFG.GT.l.) PR= (6. *X !i02/( XF.02F5* ) )**3.5*( 6./ (7.+XH02-I 
If IXFDES.CT.l. ) PRU£$»( 6 .*X«r 2 /(XMC 2 F 5 , ) )** 3,5 
* ♦{c./(7. + Xha2-:. ) )**2.5 

tC FP=Pk-OtLPR 

PFLfcS*PROci-DcLPR 
Gi TC 120 


lAdLE LOriK-OF 


7l FP=XFkl(l) 


PR INC036 
PRIN0C37 
PRIN0C38 
PRINGC3P 
PRINOGAO 
PR INOONI 
PRIN0CA2 
PRIN0CA3 
PR INCC<»A 
PR INOOA5 
PRiNOOtf 6 
))**2.5 PRINDCA7 
PRIN0LH8 
PRIN00A9 
PRIN0050 
PRIN0C51 
PRIND052 
PR IN0053 
PRINOCa'f 
PRIN0005 
PRIN001>6 
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)-XPRI( M ) / ( ) 

Jr ( J . L 1 . XKPk 1 ( * ) ) (,0 TJ 

I- 1 : : 1 = 

X! r ;T = I < n 

I- ^xrJ.G;.x^TtS^,A^.;.x'':l.l.r.x.'«PR[<I♦lJ) 
M'<.iPRI(I-fi)-XMTtRT)*(XPPI{I-fl)-<PPI(I)) 
I” ( a'’'J.'o = . X'- TE 5T . 4 'O . X iO.Lt . XfiPPI { :+l ) ) 
Cl -1 IC'j: 

r X - ' 1 ( u) + Sl p* ( X4^|J-Xl^’p I (on 

IE (>-■>•. lt. .: ) PR».: 

pA l'LS= X'-Xl { 1 ) 

1 1 (x.- D. G. L 7 . XMPRI ( 1 ) ) Gu TJ i,..; 

c i. \ n X ^ 

X ► T . 1 = X rf P P I ( I ) 


PRIN0:5 

PRInOO. 

PRInVw- 

PRINwC'- 

pK*xP* I ( I )-K xra-xxppi ( I n PR ISC' i 

PR INO^c 

GJTL'^G PRINGcc 

PRINOCb 
PR I f» :. «J C 
PR INOCo 
PR I Nouo 
PP IN0C6 
PR INO^ t, 
PKIN007 


IE ( >f- j'-i.GE .X.'irtST. ANP. Xf'OES.c ; .X'U'^I (I + i) ) PkOE!-XPKI(I ) + ( XMOLS-PR IN0c7 


/'• i - 1 ( 1 ) ) / (xrPR I ( H-. xpRi{i+*)-<pRi(i ) ) 

If (a' J S . Gi . XMTtST. and. X hJES .L c . XRPRI { I+ i n Gu TC llC 
Cl. NT n U 


PRlN0v7 
PR IN0U7 
PR1NCG7 


F-’u» G = XPPl(X-)+jLP*(AHL)Ei-XhPRI(oJ ) 
It (Pf J.- >.LT. . 1 ) PRJ :i-,l 
lie Ent:!-:. 
i^C FI T r = RR* •’RSU'J 

Ir ( 1 PR. O'! . NC . ) lKR*IPRnuM 

f T L/t' N 

ENT 


PR INir 7 
PRINjv? 
PR1NJ07 
PRIN0C7 
PRIN0C7 
PRIN0C3 
PR IKCC 6 


43 



REFERENCES 


!• Gregory, T. J.: Computerized Preliminary Design at the Early Stages of 

Vehicle Definition- NASA TM X-62,303, 1973. 

2- Ball, W- H- : Propulsion System Installation Corrections. USAF Contract 

AFFDL-TR-72-147, The Boeing Company, 1972. 

3. Sibulkin, Merwin: Theoretical and Experimental Investigation of Additive 

Drag. NACA Report 1187, 1954, 

4. Lighthill, J. J.: Higher Approximations in Aerodynamic Theory. Princeton 

University Press, 1960. 

5. Little, F. K. ; Internal Aerodynamics Manual, North American Rockwell 

Corp., Columbus, Ohio, 1968. 

6. Aerospace Handbook. General Dynamics, Fort Worth, Nov. 1962. 

7. Model Specification E1075-A, Engine, Aircraft, Turbojet. J~85-GE-13, 

General Electric Company, June 30, 1969. 

8. Crane, C. A.; Spencer, B. L.; and Vance, C. H.: Flight Manual and 

Standard Aircraft Characteristics Performance Substantiation ~ 

F-5A/B Tactical Fighter with Two J-85-GE-13 Engines. NOR-67-29, 
Northrop Corporation, Norair Division, July 1967. 

9. Kline, R. E, ; Crews, A. H. ; and Knight, W. J.: F-5A Category II Perfor- 

mance Test. Technical Report No. 65-15, Air Force Flight Test Center, 
1965. 

10. Hesse, W. J. ; and Mumford, N. V. S.: Jet Propulsion for Aerospace 

Applications. Pitman Publishing Corp., 1964. 

11. Ames Research Staff: Equations, Tables, and Charts for Compressible Flow. 

NACA Report 1135, 1953. 

12. Barry, F. W. : Conical Flow Properties for Use in the Design of Supersonic 

Inlets, Report M-1266-1, Research Department, United Aircraft Corpora- 
tion, Nov. 18, 1958. 

13. Shapiro, A. H.: The Dynamics and Thermodynamics of Compressible Fluid 

Flow, vol. I. Ronald Press Co., 1953. 

14. Anderson, W. E. ; and Nelms, Walter P., Jr.: Performance and Control of 

Air Induction Systems for a Commercial Supersonic Transport — 

1. Geometry and Steady State Performance. Presented at the AIAA/AFSC 
Vehicle Design and Propulsion Meeting, Nov. 1963. 


44 



i 



Note: A denotes areas 

(a) Inlet. 


Figure 1.- Nomenclature. 




Figure 2.- Block diagram of ACSIOIT. 
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Figure 3.- Block diagr^ of propulsion installation losses subroutine (PRINC); values computed 
are in the boxes; inputs are in parentheses belov each box. 










VO 


Installation drag coefficients: 


(D 

© 

© 


Additive = 

Bleed = 
Bypass = - 


fIl (P -P^)dAy 
^ QooAc 


■ (■'>E^E‘''Ai;{ P -P„) 

q<n Ac 


’^pV« -('n£Vg+AE(P -Peu)j 
1=0 Ac 


0 Boattail = -P»-)dAj^ 

lx. Ac 

Interference = 

1=0 Ac 


© Diverter = 

lx. Ac 

^ AatjxVcc 

(X) Auxiliary = -■ . 


Accounting method: 

Fnet = Tjjj3fp -D 

'^INST installed thrust 

D - airframe drag 
Where: 

TinST = TuniNST -(©+© < © 10+0 + ©+©)xq^ Ac 
Where: 

^UNINST uninstalled thrust (corrected for pressure 
recovery) 


Figure t-- Accounting method used in PRIWC module. 
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Figure 6,- Mass flow ratio versus Mach number for simulated F-5A with (2) 

J85-I3 engines. 
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Me.., 


(a) Maximum A/B power setting. 

Figure T-- Example propulsion installation drag calculated by PRINC for simulated F-5A w-ith (2) 

J85“13 engines; based on Pref- 






Ln 



h = 10973 m (36000 ft) 

Sref =* 15.8 (17O ft^) 

Ac = 0.083 (0.892 ft^) per engine 

Military power 


■Total interference drag 




Olxqu ‘M 


Uninstalled 



Figure 9 -- Thrust correlation for simulated. F- 5 A with (2) J85-13 engines 

h = 10973 m ( 36000 ft) . 


















Hx, 


(b) Bleed drag. 
Figure 13 .- Concluded. 
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(a) Boattail angle effects. 

Figure 15.- Correlation of nozzle boattail drag coefficient. 



Aref = 1.88 (20.2 ft^), D - I29.5 cm 

(51 in.), L = 127.0 cm (50 in.), p = 20 



f'RTTIC r.i'^dule 

Boeing data 
ref, 2 





